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ABSTRACT 

We measure the evolution in the virial mass-to-light ratio (M2oo/Lb) and virial-to-stellar mass ratio 
(M200/M*) for isolated ~ L* galaxies between z ~ 1 and z ~ by combining data from the DEEP2 
Galaxy Redshift Survey and the Sloan Digital Sky Survey. Utilizing the motions of satellite galaxies 
around isolated galaxies, we measure line-of-sight velocity dispersions and derive dark matter halo 
virial masses for these host galaxies. At both epochs the velocity dispersion of satellites correlates 
with host galaxy stellar mass, a oc M® A±0A , while the relation between satellite velocity dispersion 
and host galaxy i?-band luminosity may grow somewhat shallower from a oc L B ^ ±aA at z *~ 1 to 
a oc L° B 4±Q1 at z ~ 0. The evolution in M2Q0/-M* from z ~ 1 to z ~ displays a bimodality insofar 
as host galaxies with stellar mass below M* ~ 10 11 h^M® maintain a constant ratio (the intrinsic 
increase is constrained to a factor of 1.1±0.7) while host galaxies above M* ~ 10 11 h _1 MQ experience 
a factor of 4 ± 3 increase in their virial-to-stellar mass ratio. This result can be easily understood if 
galaxies below this stellar mass scale continue to form stars while star formation in galaxies above 
this scale is quenched and the dark matter halos of galaxies both above and below this scale grow 
in accordance with ACDM cosmological simulations. Host galaxies that are red in U — B color have 
larger satellite dispersions and hence reside on average in more massive halos than blue galaxies at 
both z ~ 1 and z ~ 0. The satellite population of host galaxies varies little between these epochs; 
the only significant difference is that satellites at z ~ 1 tend to be comparatively fainter (by ~ 0.15 
magnitudes in the mean) relative to their host luminosity than satellites at z ~ 0. The redshift 
and host galaxy stellar mass dependence of M200 /M* agrees qualitatively with the Millennium Run 
semi-analytic model of galaxy formation. 

Subject headings: galaxies: evolution — galaxies: kinematics and dynamics — galaxies: halos 



1. INTRODUCTION 

The current cosmological framework indicates that 
galaxies are embedded in massive dark matter halos that 
extend far beyond the visible baryonic component. The 
virial mass-to-light, M200/L, and virial-to-stellar mass, 
A/200 /M* , ratios provide simple distillations of the com- 
plex interplay between galaxies and their dark matter 
halos. Constraints on the dependence of M200/L and 
.M200 /M* on galaxy properties and redshift hence afford 
unique insight into the formation and evolution of galax- 
ies in a cosmological context. For example, evolution of 
M2Q0/M* provides constraints on the evolution of star 
formation efficiency in galaxies. 

Yet there are surprisingly few direct constraints on 
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M200/L or M200/A4*. Since dark matter extends far be- 
yond the visible components of a galaxy, it is in prac- 
tice difficult to probe the halo mass on large scales (i.e., 
> 100 h^ 1 kpc) due to a lack of luminous tracers. The 
halo masses of clusters of galaxies can be estimated by 
strong and weak gravitational lensing, the tight relation- 
ship between the X-ray temperature of the diffuse intra- 
cluster medium and the dark matter halo, the Sunyaev- 
Zeldovich effect, and by measuring the velocity disper- 
sion of the cluster galaxies themselves, under the as- 
sumption that the cluster gala xies are tra cing out the 
dark matter halo potential fsee lVoitll2005L for a recent 
review) . 

Probing the halos of isolated galaxies is more dif- 
ficult. The only methods that are currently ca- 
pable of directly probing the halo mass of iso- 
lated ~ L* galaxies to large radii are weak lensing 
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120051 Ivan den Bosch et alJ l2004t iConrov et all 12005b) , 
which utilizes satellite galaxies as test particles that trace 
out the dark matter halo velocity field out to several hun- 
dred kiloparsecs. One limitation of both techniques is 
that they must "stack" many isolated galaxies in order 
to build up a robust signal. Both methods have yielded 
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similar results; e.g., each finds a positive correlation be- 
tween galaxy luminosity and the virial mass of its dark 
matter halo. 

The evolution in the virial mass-to- light r atios of bright 
isolate d galaxies is only poorly constrained. Wilson ct al. 
(2001) measured the weak lensing signal for early-type 
~ L* galaxies from z = 0.8 to z = 0.1. They found lit- 
tle evolution in the halo mass of ~ L* galaxies, though 
their formal errors on the halo mass at z ~ 0.8 were 
> 50%. Furthermore, that work assumed that L* did 
not evolve with redshift and that halo mass was pro- 
portional to the square root of galaxy luminosity. Most 
recently, utilizing stellar mas ses and weak lensing data 
from COMBO-17 and GEMS. lHevmans et all l|2006|) find 
no significant evolution in the virial-to-stellar mass ra- 
tios of bright galaxies to z ~ 0.8 (though constraints are 
weak; they find that the ratio decreases by no more than 
a factor of 2.6 at la). 

Using the first ~ 25% of the recently completed 
DEEP2 redshift survey, which has now c ollected spec- 
tra for > 40, 000 galaxies at 0.7 < z < 1.4. lConrov et all 
l)2005bj) used the dynamics of satellite galaxies to mea- 
sure the halo mass for bright isolated galaxies with satel- 
lites (referred to as "host" galaxies) and compared their 
derived M/L to measurements from the Sloan Digital 
Sky Survey (SDSS) at z ~ 0. However, the small sample 
size and differences in selection effects between DEEP2 
and SDSS meant that all claims in that work had to be 
highly qualified. 

This study presents a much more detailed compari- 
son between host-satellite systems identified in the com- 
pleted DEEP2 Redshift Survey and systems found in a 
consistent way from the fourth public data release of the 
SDSS. The increased data sample at z ~ 1, combined 
with a careful handling of the different selection effects 
between the two surveys, allows a robust determination 
of the evolution in the £?-band virial mass-to-light ratio 
and the virial-to-stellar mass ratio of ~ L* host galaxies 
from z ~ 1 to z ~ 0. 

The classic Tully-Fisher relation is also capable of con- 
straining the halo mass of isolated disk-dominated galax- 
ies, though the steps required to convert the observed 
rotation speed into a dark matter halo virial mass is 
more model de pendent than either weak lensing or satel- 
lite dynamics, fvan den Boschl20 02; Kas sin et al.l f2006b; 
iGnedin et al.ll2006j) . In particular, this method must ex- 
trapolate the rotation curve far beyond the region cov- 
ered by observations and/or requires knowledge of the 
relative contribution of baryonic and dark matter to the 
observed rotation curve as a fu nction of scale (s e e, e.g ., 
Guedir, et n ,ljl2(p>. Recently. IConselice et a.U pTOffr , 
Boehm fc Zieelerl l|2006f) , and Kassin et alJ l)2006a|) have 
used the Tully-Fisher relation to constrain the evolution 
in the virial-to-stellar mass ratio from z ~ 1 to z ~ 0. 
These studies found no evidence for a change in this ratio, 
though their sample sizes were relatively small compared 
to weak lensing and satellite kinematics studies (~ 100 
objects for the first two studies and ~ 550 for the last). 

This article proceeds in the following manner: ^out- 
lines our definition for host and satellite galaxies and 
describes our method for estimating halo masses of host 
galaxies. In <JH we describe the galaxy surveys studied 
and explain the methods used to mitigate survey selec- 
tion effects, ^contains a comparison of the properties of 



host and satellite galaxies between z ~ 1 and z ~ 0, and 
in SjHlwe present satellite velocity dispersions and derived 
virial mass-to-light and virial-to-stellar mass ratios as a 
function of host galaxy redshift, luminosity, stellar mass, 
and color, ^contains a comparison between our results 
and a semi-analytic model of galaxy evolution. In Sjjjwe 
discuss these results and ^concludes. Those readers not 
interested in the technical details should focus on §§0]- 

E 

A ACDM cosmology is adopted throughout: Vl m — 
1 - fi A = 0.3, with H = 100/1 km s" 1 Mpc" 1 . All 
absol ute magnitudes quo ted here are in the AB sys- 
tem ijOke & Gunnl I1988T) . and throughout magnitudes 
are quote as Mb — 51og(/i). We adopt a mass defini- 
tion for dark matter halos such that the virial radius of 
a halo corresponds to a region with mean density 200 
times the critical density, denoted r2oo- The halo mass 
is the mass interior to r2oo, and is denoted M2oo- M* is 
reserved for stellar masses and M* for the characteristic 
scale of the luminosity function. 

2. METHODOLOGY 

This section describes our definitions of host and satel- 
lite galaxies and our methods for determining line-of- 
sight velocity dispersions and dark matter halo masses 
of host galaxies. 

2.1. Isolation and Satellite Criteria 

The use of satellite dynamics for extracting mass es- 
timates of their host galaxies is motivated by a scenario 
in which a bright galaxy resides at rest at the center of 
its dark matter halo, has no other bright companions, 
and is surrounded by faint satellites that are bound to 
the host halo and orbit within it. The criteria used to 
define host galaxies and their associated satellites are 
meant to capture such systems. Note that, according 
to our use, "host" galaxies are not completely isolated; 
rather they are isolated with respect to other compara- 
bly bright galaxies. These criteria reject objects in dense 
regions such as clusters and groups, which typically con- 
tain several comparably bright galaxies. 

Isolation criteria is specified using three parameters: 
an absolute B-band magnitude difference, AMb, a veloc- 
ity difference, AV, and a projected physical separation 
AR p . The latter two parameters define a search cylinder: 
if a galaxy has no companions within the search cylinder 
that are within AMb in absolute magnitude, then it is 
deemed to be "isolated" . 

With a set of isolated galaxies we then search for satel- 
lite companions by specifying a set of three similar pa- 
rameters: a magnitude difference, 5Mb, a velocity differ- 
ence, 5V, and a projected physical separation, 5R P (here 
and throughout we reserve A for isolation criteria and S 
for satellite criteria). We then search for companions of 
the isolated galaxies that are within the search cylinder 
defined by SV and SR P , and are fainter than the isolated 
galaxy by at least 5Mb magnitudes. Isolated galaxies 
with satellites are called "host galaxies" . 

Various authors have used different parameters for 
identifying host-satellite systems and have found that 
the resulting derived halo mass is quite insensitive to 
reasonable choices of parameters (|Prada et alJ 120031 
iConrov et alJ l2005b|) . In this work we use the search 
parameters listed in Table the same parameters are 
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TABLE 1 

Search Criteria 





Isolation Criteria 




Satellite Criteria 


Sample 


AMb 


AV 


ARp 


SM B 


SV 


SR P 


Name 




km s _1 


ft -1 kpc 




km s — 1 


h —1 kpc 


A 


1.0 


1000 


500 


1.0 


750 


350 


13 


1.5 


1000 


500 


1.5 


750 


350 


C 


1.5 


1000 


1000 


1.5 


500 


350 



Note. — The AV and SV parameters are set to large r values when 
considering samples of brighter host galaxies. See i|2,ll for details. 



used to extract systems from both the DEEP2 and SDSS 
surveys, and our fiducial set of parameters is sample A. 
In Appendix [5] wc show that our results are unchanged 
within la when adopting different search parameters. 

However, the recovered satellite velocity dispersion (see 
below) can become sensitive to SV (the maximum ve- 
locity separation between host and satellite) when the 
contribution from the true satellite dispersion is signifi- 
cant even at the edge of the dV distribution (the window 
defined by ±SV). This problem is alleviated simply by 
increasing the SV parameter for brighter host galaxies 
(which have larger satellite dispersions). Indeed, tests 
with simulations have shown that scaling the SV param- 
eter with host galaxy luminosity more robustly recov- 
ers the true satellite dispersion at large host luminosities 
Ijvan den Bosch et aLll2004l: IChen et al1l2005j) . Hence for 
the highest host galaxy luminosity and stellar mass bins 
at z ~ and z ~ 1, we increase SV by increments of 
500 km s _1 until the measured dispersion has converged 
(the satellite dispersion for fainter samples had already 
converged using sample A parameters). Convergence is 
achieved when using SV = 1000 km s _1 for all of these 
bins except for the highest stellar mass bin at z ~ 0; 
there the dispersion converges at SV = 1500 km s" 1 . 9 

2.2. Velocity Dispersion & Halo Mass Estimation 

This section describes in detail how to estimate the 
dark matter halo mass of host galaxies. This procedure 
can be conceptually separated into three steps: iden- 
tifying host-satellite systems in a galaxy catalog; re- 
constructing the line-of-sight velocity dispersion profile 
of their satellites; and converting the dispersion pro- 
file into a mass. Once host-satellite systems are found, 
the velocity dispersion profile 10 of the host galaxy dark 
matter halo can be estimated from the distribution of 
velocity differences between host and satellite galaxies, 
dV = Vhost — Ksat) m bins of projected distance from the 
host galaxy, R p . 

In an ideal world there would exist hundreds of suffi- 
ciently luminous satellite galaxies per host galaxy, and 
the task of measuring a velocity dispersion profile would 
be comparatively straightforward. Unfortunately, in 

9 In addition, for this case AV has been increased to 1500 km 
s — 1 since it is conceptually awkward for SV to be larger than AV. 

10 Note that what is actually being probed is the line-of-sight 
velocity dispersion profile; we drop 'line-of-sight' in the remainder 
of this paper for brevity. 



practice there are on average only 1 — 2 satellites per 
host galaxy, due to the magnitude limit of the redshift 
surveys used here (most isolated galaxies possess no iden- 
tifiable satellites). In order to build up a dispersion pro- 
file we must combine satellites from many host galax- 
ies and c onstruct an average profi le around an average 
host (e.g. iZaritskv k Whitelll994j) . Assuming that host 
galaxy properties such as luminosity are tightly corre- 
lated with their dark matter halo mass, then by stacking 
host galaxies in bins of lumi nosity we can reco ver the 
average underlying halo mass ( IPrada et alJEOQ^) . Weak 
lensing studies rely on the same stacking procedure, since 
the lensing signal from individua l isolated galaxies is very 
weak (e.g.. iBrainerd et alJ ll*996). 

Naively, one might expect the velocity distribution of 
satellites, f(dV), to be approximately Gaussian with a 
given by the velocity dispersion of satellite galaxie s as- 
sociated with the host galaxy ijPrada et al.1 12003^) . In 
fact, there is a significant contribution to the dV distri- 
bution from so-called "interloper" galaxies. Interlopers 
are galaxies that are classified as satellites in projection 
but are in fact not true satellites, i.e., they are not phys- 
ically associated with th e host galaxy. 

Previous authors (e.g.. iMcKav et al.l2002l IPrada et alJ 
l2003UConrov et al.ll2005b|) have modeled the effect of in- 
terlopers as a constant contribution to f(dV) at all veloc- 
ities. Tests with simulations have confirmed that mod- 
eling interlopers in this way results in a robust recovery 
of the un derlying mass distrib u tion ( see A ppendi x Q . 
Howev er, Ivan den Bosch et al.l l|2004f l and IChen et alJ 
(2005) found that in mock galaxy catalogs interlopers 
do not have a constant dV distribution but rather have 
a velocity structure somewhat similar to true satellites 
(though in the mock catalogs the width of the interloper 
distribution does not appear to scale with host galaxy 
luminosity). Despite this (or perhaps because of this), 
these authors found that modeling the interloper distri- 
bution as a constant component to the dV distribution 
does yield an accurate recovery of the velocity dispersion 
profile as well as the underlying halo mass, in agreement 
with previous work. 

Hence f(dV) is modeled as a Gaussian distribution 
plus a constant component: 



f(dV; n,a ios ) = — — + 



1 — T) 



2SV y/2itaert(5V/i/2a) 



(1) 



subject to the following normalization condition: 

+sv 

f(dV;r ] ,a loa )d(dV) = l, (2) 

sv 

where r\ is the fraction of interlopers within ±SV and o~\ os 
is the line-of-sight satellite dispersion. 

We use a maximum-likelihood method to fit this 
Gaussian-plus-constant function to the observed distri- 
bution of dV for pairs in some bin of projected separation 
R p . We maximize the likelihood function: 



hy[L( V ,o- loB )] =Y,\iky[f{dVi; 



0] 



(3) 



over a dense grid in a\ os and 77, where dVi is the dV 
value for the zth satellite-host galaxy pair, which is given 
weight Xi (see below). Since we only use dVi values for 
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pairs in some bin of projected separation, the parameters 
■q and <7i os are implicit functions of R p . 

In the following analysis satellite i is assigned a weight 
Aj according to the inverse number of satellites associated 
with the host of satellite i. This weighting scheme en- 
sures that host galaxies that have a large number of satel- 
lites (and are hence likely to be more massive than the 
average host) do not dominate the likelihood. Van den 
Bosch et al. (2004) compared weighting by host galax- 
ies (the scheme employed here) to weighting by satellite 
galaxies (which would be equivalent to setting A = 1 for 
all satellites) in simulations and found that the recovered 
velocity dispersion differs between these two schemes as a 
function of host luminosity 11 . Weighting by host galaxy 
more fairly represents the average mass of host galaxies 
within a given luminosity or stellar mass bin. 

The resulting fit yields not only a measurement of the 
velocity dispersion but also an estimate of the interloper 
fraction. Marginalizing the likelihood over rj provides an 
estimate of the ler errors on the velocity dispersion. 

Redshift uncertainties are accounted for by subtract- 
ing in quadrature the rms redshift error, tr e rr, from the 
measured velocity dispersion, a\ os : 

°cst = ^<jf os -2al,.. (4) 

The redshift uncertainty enters twice because we are sub- 
tracting two velocities, the host from the satellite. The 
resulting velocity dispersion, er ost , is then our best esti- 
mate of the true dispersion. For the range of velocity 
dispersions probed here, folding redshift errors from the 
DEEP2 and SDSS surveys (cr crr < 30 km s" 1 ) into our 
measured dispersion changes results by only a few per- 
cent. In the following sections we simplify our notation 
to a = cr ost for brevity. 

With an estimate of the velocity dispersion of the 
host galaxy dark matter halo we can, with addi- 
tional assumptions, extract th e virial mass of the halo. 
We follow the proc edure of IPrada et alJ l|2003|) and 
IConrov et alJ l|2005b(l in converting velocity dispersions 
to virial masses. 

The density profile of a dark matter hal o is parameter- 
ized using the NFW ijNavarro et al.llT997j) model: 

pi (r/r s )(l + r/r s f [> 

where p° c is the present critical density, r s = r2oo /c, and 

200 c 3 
3 ln(l + c) - c/(l + c) ' [ } 

where r2oo is defined as the radius where the mean in- 
terior density is 200 times the critical density. The 
concentration, c, is inversely related to the mass of a 
dark matter halo, and, at fixed mass, scales with red- 
shift as (1 + z)- 1 ijBullock et alJl200H) ; concentrations 
for dark matter halos hosting galaxies or clusters range 
from 3 < c < 25. In our analysis we fix c = 10, which 
is consistent with the concentration of a ~ 10 12 /i _1 M 
halo at z = 0, and scale c by (1 + z)^ 1 for higher redshift 

11 This is due to the fact that brighter galaxies will have more 
satellites, especially in a flux limited sample, and hence the dif- 
ference in weighting schemes is more pronounced in the regime of 
bright host luminosity. 



samples. However, as demonstrated below, the density 
profile depends only weakly on concentration over the 
scales probed (20 < Rp < 150 h^ 1 kpc), and hence the 
assumed concentration has little effect on the resulting 
mass estimates (cf. Appendix lB|l . Throughout we as- 
sume that the satellite galaxies follow the radial profile 
of the dark matter. The effect of a spatial bias between 
satellites and dark matter on the recover ed mass have 
been show by Ivan d en Bosch et al. ( 2004) to be at the 
few percent level (see also Appendix . 

We then determine the radial velocity dispersion pro- 
file by integrating the Jeans equation, which relates 
the density profile and gravitational potential to the 
radial velocity dispersion, using Equation [S] to specify 
the potential. Finally, we integrate along the line-of- 
sight. Both integrations require knowledge of the veloc- 
ity anisotropy, = 1 — of/o-j_, of the satellite popu- 
lation. Fortunately the line-of-sig ht velocity dispersion 
profile depends only weakly on l)van den Bosch et alJ 
l2004UMamon fc LokarJElO-Vi. In AppendixEwe explore 
both an isotropic distribution ( = 0) and an ani s otrop y 
parameterization suggested bv iMamon &: Lokasl {2005), 
and confirm that the derived mass is robust to assump- 
tions about 0. Below we set = 0. 

Given these assumptions, there remains only one free 
parameter, the normalization of the dispersion profile, 
which is related to the mass within r200j denoted A/200 ; 
of the dark matter halo. The normalization is obtained 
via x 2 minimization using the measured a(R p ) points. 
The majority of the analysis below uses only one velocity 
dispersion measurement to derive a virial mass. In Ap- 
pendix [B] we demonstrate that estimating virial masses 
with only one dispersion measurement does not bias the 
recovered mass. Including more velocity dispersion mea- 
surements in finer projected separation bins simply has 
the effect of decreasing the error on the recovered mass (if 
there are sufficient numbers of satellites to increase the 
nu mber of radial bins, w hich is not the case in DEEP2). 
See iKlvoin et al.1 (|2006) for a detailed study of the full 
satellite velocity dispersion profile measured for SDSS 
host galaxies at z ~ 0. 

3. THE DATA 

This section presents the low- and high-redshift galaxy 
catalogs used to identify host and satellite galaxies, and 
describes how to account for the differences in selection 
effects between the two catalogs. 

3.1. The SDSS 

The Sloan Digital Sky Survey (SDSS:lYork et al.ll2000t 
lAbazaiian et all20()llAdelman-McCarthv et al.l2006[) is 

an extensive photometric and spectroscopic survey of the 
local Universe. As of Data Release 4, imaging data ex- 
ist over 6670 deg 2 in five bandpasses, u, g, r, i, and 
z. Approximately 670, 000 objects over 4780 deg 2 have 
have been targeted for follow-up spectroscopy as part 
of SDSS are included in DR4; mo st spectroscopic tar - 
gets are brighter than r = 17.77 ijStrauss et alj [2002). 
Automated software performs all the necessary data re- 
duction including the assignment of redshifts. Redshift 
errors are < 30 km s^ 1 , similar to DEEP2. The spec- 
trogra ph tiling algorithm e nsures nearly complete sam- 
pling l|Blanton et alJl2003aD . yet the survey is not 100% 
complete due to several effects: 1) fiber collisions do not 
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allow objects separated by < 1' to be simultaneously tar- 
geted, resulting in ~ 6% of targetable objects failing to 
be targeted for spectroscopy; 2) a small fraction (< 1%) 
of targeted galaxies fail to yield a reliable redshift; and 
3) bright Galactic stars block small regions of the sky. 
None of these effects is expected to impact our analysis. 
The overall completeness of the SDSS, as defined by the 
number of objects with successful redshifts divided by the 
number of objects in the imaging catalog with r < 17.77, 
is ~ 90%. The parent catalog has 166, 923 high quality 
redshifts between 0.01 < z < 0.1 and 100° < RA < 275°. 

For this analysis we make use of the hybrid N YU Value 
Added Galaxy Catalog (VAGC) lIBlanton et a,ll2005D . In 
addition we use the publicly available package kcorrect 
v4.1.4 l|Blanton et al.l l2003ht iBlanton fc Roweisl 12006(1 
to derive restframe i?-band magnitudes and U — B colors 
for SDSS galaxies. All SDSS galaxies are K-corrected to 
z = 0.0. Galaxies are divided into red and blue popula- 
tions based on the valley visible i n the color-magnitude 
diagram (e.g.. iBaldrv et aTll2004l) . We account for the 
fact that the valley moves redward for brighter galaxies 
with the following color-cut: 



U-B = -0.066 M B -0.05. 



(7) 



We also obtain stellar masses for SDSS galaxies using 
kcorrect v4.1.4 routines. These stellar masses, which 
have been obtained assuming a Chabrier IMF, are con- 
sistent with the stellar masses of lKauffmann et all |2003) 
but are lower by ~ 0.3 dex f r om th e color-based stellar 
mass estimates of iBell et all (12003 i. The offset is due 
primarily to differences in the assumed IMF. 

3.2. The DEEP2 Survey 

Th e DEEP2 Galaxy Redshift Survey l|Davis et all 
2003J has gathered optical spectra for ~ 40, 000 galax- 
ies in the redshift ran ge 0.7 < z < 1.4 using the 
DEIMOS spectrograph l|Faber et ail 120031) on the Keck 
II 10-m telescope. The survey spans a comoving vol- 
ume of ~ 5 x 10 6 h~ 3 Mpc 3 , covering ~ 3 deg 2 over 
four widely separated fields. Target galaxies were se- 
lected using BRI imaging from the CFHT telescope 
down to a limiting mag nitude of R = 24.1 ijCoil et al.l 
2004). In three of the four fields we also use apparent 
colors to exclude objects likely to have z < 0.7. This 
pre-selection greatly enhanc es our efficiency fo r target- 
ing galaxies at high redshift l|Faber et al.ll200^) . Due to 
the high spectral resolution (R ~ 5, 000) and excellent 
sky subtraction provided by the DEIMOS spectrograph 
and DEEP2 /DEIMOS data reduction pipeline (Cooper 
et al., in prep), our rms redshift errors are ~ 30 km s _1 
determined from repeated observations. Details of the 
DEEP2 observations, cat alog constru c tion, a nd data re- 
duction can be found in IDavis et all l)2003|) . iCoil et all 
(2004), and IDavis etafl l|2005jb Restframe U — B colors 
and absolute .B-band magnitudes have been derived as 
described in IWillmer et all (|2006f) . The parent DEEP2 
catalog includes 21, 184 galaxies with high-quality red- 
shifts in the redshift interval 0.70 < z < 1.2. 

Stellar masses for DEEP 2 galaxies are estimated in 
the following way. For the subset of DE EP2 galaxies 
for whi ch there exists if s -band imaging, IBundv et alJ 
( 2005b) has determined stellar masses b a sed o n the 
methodology outlined in iKauffmann et alJ l)2003|) with 
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Fig. 1. — Color-magnitude diagram for host galaxies (blue/black 
points) and satellites (green/grey points). Galaxies drawn from 
the SDSS at z ~ are plotted in the top panel, while the bottom 
panel shows galaxies at Z ~ 1 from the DEEP2 survey. The top 
panel contains only 25% of the total number of objects for clarity. 
The dashed broken line defines the completeness limit at z = 1 
(bottom panel), and a similar limit at z = (top panel), where 
in this case the line has been shifted to the right according to the 
estimated evolution in Mg from z = 1 to z = 0. The dotted line 
in each panel indicates our division between red and blue galaxies. 



a Chabrier IMF. We then use an empirically derived re- 
lation between restframe UBV colors and stellar mass 
(C.N. A. Willmer, private communication) to obtain stel- 
lar masses for DEEP2 galaxies that do not have if s -band 
imaging. The stellar masses obtained in this way agree 
well with the stellar mass estimates obtained for DEEP2 
galaxies from the kcorrect v4 . 1 . 4 routine (cf. the pre- 
vious section), as expected since both methods use the 
same IMF. Since the same IMF is used, we do not expect 
any systematic offset between the stellar mass estimates 
at z ~ and z ~ 1. 

As with the SDSS data, DEEP2 galaxies are split into 
red and blue populations based on the valley in the color- 
magnitude diagram. Th ese two popula t ions a re divided 
in a manner identical to IWillmer et al.l <)2006j) using the 
following color-cut: 



U — B = -0.032 (M B + 21.63) + 1.03. 



(8) 



The DEEP2 survey spectroscopically targets ~60% of 
objects that pass the apparent magnitude and color cuts 
mentioned above. Of those targeted objects, we are able 
to secure redshifts for > 70%. Follow-up observations 
have shown that ~ 15% of the targets are objects at z > 
1.5 and fail to yield redshifts from DEEP2 for that reason 
(C. Steidel, private communication). We therefore have 
successful redshifts for 60%x85%= 51% of all galaxies 



6 



% 0.00 



. Overall 

■ Hosts 

'- 




»- 


y 

4 

i 

■ 


z~0 ~ 


i 


'i''' 
; - 

>m \ 


• i 

i 


z~l . 



0.0 0.5 1.0 1.5 2.0 

U-B 

Fig. 2. — U — B colors for host galaxies (solid lines) compared 
to all galaxies of similar luminosities and redshifts (dashed lines), 
for galaxies at z ~ (top panel) and z ~ 1 (bottom panel). The 
comparison sample was restricted to have the same rcdshift and 
luminosity distributions as the host galaxy sample. The histograms 
are plotted in units such that the integral under each curve is one. 



at 0.7 < z < 1.4 in the surveyed fields with apparent 
magnitude of R < 24.1. 

3.3. Defining Uniform Samples from SDSS and DEEP2 

Every galaxy survey is unique; unfortunately, that 
makes comparison between surveys difficult. In our case, 
there are two separate effects that must be taken into ac- 
count in order to compare the DEEP2 and SDSS galaxy 
surveys fairly. The first issue is the differing sampling 
rates; this is trivial to resolve, however. As mentioned 
in previous sections, the completeness of the DEEP2 
and SDSS surveys are ~ 50% and ~ 90%, respectively. 
This difference, if not accounted for, would result in a 
much larger fraction of falsely isolated host galaxies in 
the DEEP2 sample (i.e., galaxies that appear isolated 
in the spectroscopic catalog but are not truly isolated) 
and could bias the resulting satellite velocity dispersion 
profile with respect to the SDSS sample. 

In order to mitigate this difference we randomly di- 
lute the SDSS redshift catalog to the same complete- 
ness as DEEP2 (that is, 40% of SDSS galaxies are ran- 
domly removed from the catalog). Extensive tests with 
mock catalogs have confirmed that for the DEEP2 survey 
the incompleteness is very close to uni form as a function 
of 3-dimensional ga laxy over-density ijGerke et alJl2005l 
ICooper et al.ll2005|) . hence this simple random dilution 
is sufficient for our purposes. As shown in Appendix lAl 
this dilution results in a ~ 30% increase in the estimated 
virial masses of host galaxies compared to the complete 



(undiluted) SDSS. This should be kept in mind when 
considering the results in 

The second effect primarily impacts the satellite pop- 
ulation. While both surveys select targets that are 
brighter than an apparent i?-band magnitude limit, in 
the DEEP2 survey this selection corresponds to an ap- 
proximately restframe _B-band selection at z ~ 0.7 and 
approximately restframe {7-band at z ~ 1, while in the 
SDSS this selection corresponds closely to restframe R- 
band. This means that, for the DEEP2 survey, as one 
considers fainter objects, redder galaxies will drop out of 
the survey at a brighter Mb than bluer galaxies. Thi s 
selection effect is well understood ijWillmer et alj f20061. 
and is accounted for it in the following way. 

Figure ^ shows the color-magnitude distribution for 
all host galaxies (blue/black points) and their associ- 
ated satellites (red points). Results for host and satellite 
galaxies at z ~ are shown in the top panel, while re- 
sults for z ~ 1 are in the bottom panel. The dotted lines 
indicate our division into red and blue populations. Note 
the increase in the number of red host galaxies between 
z ~ 1 and z ~ at bright luminosities, where our data 
are complete. The increasing abundance of r ed galaxies 
with time has been studied in detail elsewher e ijBell et alJ 
l200llFaber eTZll2(M IWillmer et alJ l2006). 

At a specified redshift, the apparent DEEP2 i?-band 
magnitude limit can be modeled by a broken line in the 
co lor-magnitude diagra m (cf. Fi gure U here a nd Fig ure 4 
in lWillmer et alJ <|200fift ; see also Gcr ke et all l|200fih for a 
more detailed discussion of this effect). In particular, the 
broken line in Figure ^ can be used to define a volume 
limited sample at z < 1 that uniformly samples color- 
magnitude space; when selecting such samples we use 
only DEEP2 galaxies with z < 1 and include only galax- 
ies that are brighter than this line in color-magnitude 
space. Note that the broken line is a function of rcdshift. 
Moreover, the dearth of faint red galaxies in the bottom 
panel is attributable to this i?-band selection effect. For 
the sample at z ~ the broken line is shifted accord- 
ing to the observed evolution in the B-band luminosity 
function, M B (z) = M B (z = 0) — 1.37z, which is approx- 
imately independent of color ijFaber et alJl2005|) . Thus, 
selecting galaxies brighter than this broken line ensures 
that SDSS and DEEP2 are similarly complete relative 
to M B as a function of both luminosity and color at all 
redshifts z < 1. 

For host galaxies this selection effect is not particularly 
important because host galaxies are all rather luminous 
(by definition, since these galaxies must be bright enough 
to have satellite galaxies that are at least 8Mb magni- 
tudes fainter). For such bright galaxies both very red 
and very blue host galaxies are included in both surveys 
over the full rcdshift ranges we consider. This is appar- 
ent in Figure Q where nearly all of the blue/black points 
in both panels (representing host galaxies) are to the left 
of the broken dashed line. 

Unfortunately, accounting for the apparent i?-band se- 
lection effect drastically reduces the number of avail- 
able satellite galaxies, as can been seen in Figure ^ 
where most of the green/grey points (representing satel- 
lite galaxies) are to the right of the broken dashed line. In 
order to increase our statistics, when measuring velocity 
dispersions in |S] we do not account for this selection ef- 
fect in the manner mentioned above. Assuming that the 
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Fig. 3. — Fraction of satellites as a function of dM = M|j at — 
Mhost in sample A. Satellite-host systems at z ~ (solid line) 
are compared to systems at z ~ 1 (dashed line). Selection differ- 
ences between the parent catalogs from which these systems were 
extra cted are carefully accounted for using the "i?-band cut" (see 
for details). 



measured velocity dispersion does not depend on satellite 
properties, then the velocity dispersion should be insen- 
sitive to this selection effect, since including this effect 
will only decrease the total number of host galaxies but 
will not preferentially select one type of host galaxy (e.g., 
bright or red) over another 12 . However, when comparing 
host and satellite properties between DEEP2 and SDSS 
in 21 we account for this selection effect since it strongly 
impacts the satellite population. 

4. PROPERTIES OF HOST GALAXIES AND THEIR 
ASSOCIATED SATELLITES 

In this section we present several salient properties of 
host and satellite 13 galaxies and investigate the evolution 
of these properties from z ~ 1 to z ~ 0. Differences in 
the selection effects in DEEP2 and SDSS are taken into 
account as described in the previous section. 

Table lists the different search criteria used for iden- 
tifying satellites and host galaxies. Our fiducial sample is 
A; this section and the next contains results using those 
search criteria only. These search p arameters are simi- 
lar to ones adopted in previous work ijMcKav et alJl2002T: 

12 The situation may not be this simple if host galaxy proper- 
ties are strongly correlated with satellite properties. For example, 
if red satellites preferentially exist around red host galaxies, then 
by missing red satellites we would be implicitly missing red host 
galaxies. Such a correlation between central and satellite galaxy 
prope rties has recently been observed at z ~ iWcinmann et al. 
2006), though the signal is not large for the types of systems ex- 
plored here. Indeed, the fraction of red host galaxies does not 
significantly change when including or neglecting the ii-band ef- 
fect. 

13 Here we are actually presenting the ensemble properties of 
true satellites and interloper galaxies. Unless explicitly stated oth- 
erwise, we refer to the combination of true satellites and interlopers 
as "satellites" hereafter. 



iPrada et al.ll200l iConrov et al.ll2005b|) . In Appendix |U 
we demonstrate that our results are robust to the partic- 
ular choice of parameters. 

Table |2 contains information on the host and satellite 
samples using search criteria A. The Table includes the 
redshift and magnitude ranges over which we search for 
hosts and satellites, the median redshift of each sample 
of host galaxies, the number of hosts and satellites found 
in each sample, and the mean satellite luminosity. This 
information is tabulated for samples restricted to the re- 
gion of color-magnitude space where the parent catalogs 
are complete at both epochs (the samples used in this 
section, labeled in the table "with i?-band cut"), and 
samples that include regions of color-magnitude space 
where the parent catalogs are not complete (the samples 
used in the following section, labeled "without i?-band 
cut" in the table). 

In order to quantify the distribution of U — B colors for 
host galaxies, Figure [3 plots histograms aiU — B color at 
z ~ 1 (top panel) and z ~ (bottom panel) . Each distri- 
bution is normalized such that its integral is unity. The 
sample of host galaxies (solid lines) is compared to an 
"overall" sample of galaxies (dashed lines) . The compar- 
ison sample has been constructed such that it samples 
the redshift-luminosity plane with the same density as 
the host galaxies. It is evident that host galaxies have a 
distribution in U — B colors comparable to all galaxies 
at the same luminosity at both z ~ 1 and z ~ 0, and 
that host galaxies at z ~ are redder than host galaxies 
at z ~ 1 only insofar as the overall galaxy population 
is redder at z ~ compared to z ~ 1. One can also 
clearly see the growth in the abundance of red galaxies 
between z ~ 1 and z ~ 0, as noticed in previous studie s 
(|Bell et alJl200l iFaber et al]l2"F)05t iWillmer et all2f)f)^ . 

The samples used in Figure |21 have the i?-band cut 
taken into account. However, the full samples (neglecting 
the i?-band cut) display nearly identical host U ~B color 
distributions when compared to the i?-band cut samples. 
This indicates that the full samples, which are used in SJ5J 
do not contain any artificial changes in the U — B color 
distribution of galaxies with time. 

Generating a fair comparison sample is a requisite for 
interpreting the distribution of colors of host galaxies. 
Since host galaxies are in general much brighter than an 
average galaxy (i.e., a galaxy drawn at random from the 
full survey), without a fair comparison sample we would 
have falsely concluded that host galaxies are redder than 
the average galaxy. These selection effects are exacer- 
bated in the DEEP2 sample, where the red galaxy pop- 
ulation is a strong function of redshift due both to the 
apparent i?-band limit of the survey and to evolution in 
red number density, but are effectively mitigated with a 
proper comparison sample. 

Figure El plots the distribution, f s (dM), of magnitude 
differences, dM = Mff* - M| ost at z ~ 1 (dashed line) 
and z ~ (solid line). Error bars denote la Poisson 
uncertainties. We conclude from these distributions that 
the average host galaxy at z ~ 1 has fainter satellites 
than an average host galaxy at z ~ 0. This conclusion 
is unchanged if we separately consider f s (dM) in bins of 
host galaxy luminosity. The average dM for each sample 
reflects this difference as well, though less strikingly: at 
z ~ 1 (dM) = 1.60 while at z ~ (dM) = 1.46. 



Note that since these distributions of f s (dM) are nor- 
malized to the total number of satellites at each redshift, 
these results are insensitive to the presence of interlop- 
ers unless interlopers have a dM distribution that varies 
with redshift. Although such a scenario seems particu- 
larly nefarious, it cannot explicitly be ruled out. 

In addition we have investigated the sensitivity of these 
results to our assumed evolution in the luminosity func- 
tion. Throughout we have assumed that the characteris- 
tic scale of the luminosity function evolves as M%(z) = 
M% (z = 0) - 1.37 z independent of color ijFaber et all 
2005). The evolution in the luminosity function is im- 
portant here because it determines how the dashed line 
in Figure ^ evolves with redshift, which in turn defines 
the samples used in this section. If the amount of evo- 
lution in Mi is varied by ±0.3z (which is t he la uncer- 
tainty in the evolution of M B as reported bv lFaber et al.1 
2005) , the qualitative result of Figure namely that 
DEEP2 satellites are on average fainter with respect to 
their host luminosities than SDSS satellites, remains un- 
changed. However, the case of maximal evolution, where 
M B evolves by 1.67 magnitudes per unit redshift, results 
in a significantly less striking difference at large dM, for 
the obvious reason that this maximal evolution in M B 
allows many more faint satellite galaxies to be included 
in the sample at z ~ (in essence, the dashed line in the 
top panel of Figure moves to the right). Thus, if the 
difference in satellite properties seen in Figure [3] is not 
re al, then evolution i n M B is even stronger than found 
bv lFaber et"aD {20051 . 

Moreover, dynamical friction acting on the z ~ 1 host- 
satellite population would tend to produce a trend oppo- 
site to what is observed here. Because dynamical friction 
is most efficient when it acts between objects of compa- 
rable mass, it will cause the brightest satellites to sink 
toward the host preferentially, resulting in comparably 
more fainter satellites at z ~ than at z *~ 1, contrary 
to our observations. 

In Figure 0] we plot the fraction of host galaxies with 
at least iV sat satellites, f s (> N sat ), as a function of N sat , 
both at z ~ 1 (dashed line) and z ~ (solid line). Er- 
rors reflect Poisson uncertainties. It is evident that, when 
comparing volume limited samples at z *~ 1 and z ~ 
which are similarly complete in color-magnitude space, 
host galaxies at high redshift are associated with more 
satellites than host galaxies at low redshift (although the 
difference is weak, see below). This is illustrated alter- 
natively by considering the average number of satellites 
per host galaxy: (N) — 1.26 at z ~ 1 and (N) = 1.15 at 
z~ 0. 

In this case interlopers have a more direct impact. In 
jJSl we show that the interloper fraction decreases from 
~ 21 ± 6% at z ~ 1 to ~ 16 ± 4% at z ~ (so that 
within la the interloper fraction is constant with red- 
shift). Therefore, the changes in the average number of 
satellites per galaxy could be due to changes in interloper 
contamination with z at < la. Uncertainties in the evo- 
lution of M B from z ~ 1 to z ~ also strongly impact 
the results shown in Figure 01 because a larger/smaller 
evolution in Mb than what is assumed here will result 
in more/less faint galaxies in the z ~ "i?-band cut" 
sample, which will in turn result in more/less satellites 
at z r-j compared to z ~ 1. In short, these results are 
unfortunately too sensitive to various assumptions and 
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Fig. 4. — Fraction of host galaxies that have > -/V sa t satellites, as 
a function of AT sa t ■ As in the previous figure, we use search criteria 
A and account for differences between the parent c atalo gs at z ~ 1 
and z ~ using the "i?-band cut" as described in H3.3I 



TABLE 2 

Summary of Samples 





SDSS 


DEEP2 


redshift range 


0.01 < z < 0.10 


0.7 < z < 1.2 


<*) 


0.06 


0.84 


Host Mg — 51og(/i) range 


< -19 


< -20 


without R-band cut: 






Total Sample Size 


102, 656 


21, 184 


Satellite {M B - 51og(/i)) 


-18.3 


-19.4 


iV" sa t 


5414 


1007 


-^host 


3431 


755 


with R-band cut: 






Total Sample Size 


56, 397 


12,887 


Satellite (M B - 51og(/i)) 


-19.0 


-19.6 




1475 


554 




1283 


440 



Note. — The host and satellite galaxies contained in these 
samples were obtained using search criteria A (cf. Table HI . 
The total sample size for the SDSS survey refers to the survey 
after it has been diluted by 40% to match the completeness 
of DEEP2. 



uncertainties to provide robust conclusions regarding the 
differential evolution in satellite numbers between z *~ 1 
and z ~ 0. 

5. VELOCITY DISPERSION MEASUREMENTS & DERIVED 
VIRIAL MASSES 

We now present the measured satellite velocity disper- 
sion, a, as a function of host galaxy redshift, luminos- 
ity, stellar mass, and color. In addition, we derive virial 
mass-to-light ratios (M2oo/Lb) and virial-to-stellar mass 
ratios (M200 / 'M m ) at z ~ 1 and z ~ 0. In this section the 
i?-band selection effect is not accounted for when com- 
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paring samples at z ~ 1 and z ~ 0; see i j3.3l for details. 
Results are for host and satellite galaxies identified ac- 
cording to search criteria A (cf . Table . Since both 
SDSS and DEEP2 are constructed to have the same com- 
pleteness (~ 50%), all virial masses quoted herein are 
~ 30% higher than they would be had 100% complete 
surveys been used. This is attributable to systems that 
only appear isolated in the diluted sample, but in fact 
reside within more massive halos with multiple bright 
companions (see Appendix . 

Appendix 151 demonstrates that the results presented in 
this section are robust to various assumptions, including 
the specific search criteria used to define the samples and 
the velocity anisotropy and concentration of the satellites 
orbiting within their host galaxy's dark matter halo. Ap- 
pendix O presents an assessment of and motivation for 
the broader methodological assumptions inherent in us- 
ing the motions of satellite galaxies to extract host galaxy 
halo masses. 

The nominal redshift limit of the SDSS parent catalog 
has been extended from z = 0.1 to z = 0.2 for the highest 
stellar mass bin, doubling the number of satellites in that 
bin. Separately measuring the dispersion in this stellar 
mass bin for the fiducial sample with z < 0.1 and the ex- 
tended sample with z < 0.2 results in differences within 
10 km s . In other words, adding these higher redshift 
data does not appear to bias the resulting measurement, 
but it significantly decreases the errors on the dispersion 
due to the increased number of satellites. Increasing the 
redshift limit of the parent catalog for the other sam- 
ples has little effect since the apparent magnitude limit 
of the SDSS yields few faint satellite galaxies at higher 
redshifts. 

5.1. Results as a Function of Host Galaxy Luminosity 
and Stellar Mass 

In Figurc[S]we present the relative velocity of satellites, 
dV, as a function of host galaxy stellar mass, M*, for all 
satellites within R p = [20, 150] h^ 1 kpc at z ~ (top left 
panel) and z ~ 1 (top right panel) 14 . Note again that 
the samples used in this and the following figures do not 
account for the different survey selection effects between 
DEEP2 and SDSS (which we argued in JOD should not 
impact the results in this section) and is hence a super- 
set of the samples used in 21 O ne can see clearly that the 
satellite velocity dispersion is increasing with increasing 
host stellar mass. The solid vertical lines indicate the full 
width at half maximum given by our dispersion measure- 
ments, in bins of host stellar mass. The bottom cluster 
of panels shows the distribution of relative satellite ve- 
locities in bins of stellar mass. The smooth lines indicate 
our best fit Gaussian+constant models for host weight- 
ing (solid lines) and satellite weighting (dotted lines). 
The difference between weighting schemes becomes in- 
creasingly important for higher stellar mass bins, since 
the hosts in these bins have on average more satellites. 
We only display results as a function of stellar mass for 

14 As we show in Appendix |E] the results presented in this sec- 
tion do not change significantly when using larger, smaller, or more 
bins in R p . Our fiducial bin size is motivated by the fact that 
interloper fractions increase with increasing R p , so although one 
includes more true satellites with a larger maximum R p , the larger 
interloper fraction results in no significant improvement in the dis- 
persion measurement. 



TABLE 3 

Velocity dispersions and halo masses of host galaxies as 

A FUNCTION OF Mb 



Ma bin 



(M B ) color N Bl 



<y M 2 oo/Lb 

km s _1 HMq/Lq 



[-19.25, 


-19.75] 


-19.5 


[-19.75, 


-20.25] 


-20.0 


[-20.25, 


-20.75] 


-20.5 


[-20.75, 


-21.50] 


-19.5 


[-19.25, 


-20.00] 


-19.7 


[-20.00, 


-20.50] 


-20.3 


[-20.50, 


-21.50] 


-20.9 


[-19.25, 


-20.00] 


-19.7 


[-20.00, 


-20.50] 


-20.3 


[-20.50, 


-21.50] 


-20.8 



Results at 
all 
all 
all 
all 
red 
red 
red 
blue 
blue 
blue 

Results at z 



-19.50, 


-20.75] - 


20.4 


all 


-20.75, 


-21.50] - 


19.5 


all 


-21.50, 


-23.00] - 


21.9 


all 


-20.50, 


-22.00] - 


21.3 


red 


-20.50, 


-22.00] - 


21.2 


blue 


Note. 


— Absolut 


e B-band ma 



~ 0: 

332 

505 

657 

413 

299 

366 

475 

267 

247 

215 

~ 1: 

79 

154 

117 

133 

149 



132 


Ml 
-ii 


228!^ 


144 


Ml 
-9 


1801™ 


188 


M3 
-11 


254lj* 


235 


h-6 
-24 


277lii* 


15913 


00 '-150 


215 


M5 
-IS 


" t0 ' i -141 


2581^ 


o 84 +143 
J °^-124 


118 


f 11 
-11 


166+ff 


139 


M3 

-11 




186 


M8 
-18 


188l D 8 



H8lf 8 
1531^ 
272129 
23ll« 
1441?° 



--,+45 
' — 37 

63±i* 

13011? 

15ll 7 « 

* 3 -21 



ides, Mb, are quoted as 
Mb — 51og(/i). Recall that masses and dispersions are systemat- 
ically overestimated (by ~ 30% in mass) due to incompleteness 



effects (cf. Appendix and i|3.3l 



TABLE 4 

Velocity dispersions and halo masses of host 
galaxies as a function of m, 



M* bin (M,) N Si 



km s 



M 2m /M t h 



[9.8, 10.4] 
[10.4,10.7] 
[10.7,11.0] 
[11.0,11.5] 

[9.6, 10.4] 

[10.4,11.0] 

[11.0,11.6] 



Results at z ~ 0: 



10.2 
10.6 
10.8 
11.1 



422 
659 
663 
483 



H8l 9 9 
1551 11 
197115 
423l| 



105l 2 31 
83lH 
8llf 9 
333113 



Results at z ~ 1: 

10.1 74 109l?g 
10.7 153 162l?n r >7 

11.2 121 290 



,+24 
-20 
-,+31 
-40 



+61 
-35 
+ 2!) 
2X 

77+46 
"-30 



Note. — All stellar masses, M*, are quoted in 
units of log(M* h~ 2 Mq). Recall that masses and dis- 
persions are systematically overestimated (by ~ 30% 
in mass ) du e to incompleteness effects (cf. Appendix 
|A]and JOl . 



simplicity, but we also fit for a in bins of Mb- 

Figure contains our main results. The top panels 
present the satellite velocity dispersion measured within 
a projected separation of R p = [20, 150] hr 1 kpc as a 
function of £>-band magnitude (top left panel) and stel- 
lar mass (top right panel) for galaxies at z ~ 1 (blue 
diamonds) and z ~ (black circles). Thee la errors on 
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Fig. 5. — Left panels: The top panel plots the relative line-of-sight satellite velocity, dV, as a function of host stellar mass, M*. for all 
satellites within R p = [20, 150] h~ 1 kpc at z ~ from sample A. The solid vertical lines indicate the full width at half maximum taken 
from the measured dispersion within four bins in host stellar mass. The bottom panels display the dV distribution in bins of stellar mass 
normalized to the total number of satellites within the bin (stellar mass bin width is displayed in the top left corner in units of log(M» 
h~ 2 Mq)). The top right corner notes the number of satellites within each bin. We over-plot the best fit Gaussian+constant model weighted 
by host galaxy (used in the following analysis; solid lines) and weighted by satellite galaxy (dotted lines). Right panels: Same as the left 
panels, now for data at Z ~ 1. 



the dispersion measurement were obtained from the like- 
lihood contours described in and error bars in the 
B-band magnitude on the x-axis represent the 68% range 
within each magnitude bin. The information in these fig- 
ures is also listed in Tables|21and0]where we list, for each 
bin in host Mb and M*, the number of satellites, mean 
Mb and M* of host galaxies, recovered satellite velocity 
dispersions, and virial-to-stellar mass and mass-to-light 
ratios. 

In the bottom panels we present the same informa- 
tion, where now the velocity dispersion measurements 
have been converted into virial masses (cf. ij2.2J) and 
quote our results in terms of the virial mass-to-light ra- 
tio (M200/LB, bottom left panel) and the virial-to-stellar 
mass ratio (M 2 oo /M* , bottom right panel) . Note that in 
order to convert the velocity dispersion into a virial mass 
we use only the velocity dispersion within R p — [20, 150] 
h^ 1 kpc rather than attempting to measure the disper- 
sion in multiple bins; this is due to our limited statistics, 
especially at z *~ 1. 

Errors on the virial mass-to-light ratios reflect the er- 
rors on the virial mass, which were obtained from \ 2 
minimization between the observed velocity dispersion 
and the theoretical dispersion profile. Any errors in the 
luminosity are not included when calculating the error 
on M200/LB] we simply use the mean luminosity within 



each bin (the median £?-band magnitude agrees with the 
mean to within ±0.01). Error bars on the B-band mag- 
nitudes on the x-axis again represent the 68% range in 
each magnitude bin. The virial-to-stellar mass ratios are 
plotted similarly. 

We mention in passing that the mean redshift of host 
galaxies does not vary strongly across the B-band mag- 
nitude and stellar mass bins. At z ~ the mean redshift 
varies from 0.044 in the faintest bin to 0.076 in the bright- 
est bin, while at z ~ 1 it increases from 0.80 to 0.91 be- 
tween the faintest to brightest magnitude bins explored 
here. 

In addition to velocity dispersion measurements, our 
maximum-likelihood technique provides an estimate of 
the interloper fraction (the number of interlopers divided 
by the total number of satellites) for each sample. The 
interloper fraction for the z ~ samples is ~ 16 ± 4% 
and is constant (within la) across the host luminos- 
ity and stellar mass bins. This interloper fraction at 
z ~ is in agreem ent with previous work. In particular, 
IPradaeTaO J2003) found interloper fractions of 17-20% 
depending on their sample definition. At z *~ 1 the inter- 
loper fractions are noisier due to the smaller number of 
satellites. We measure an average interloper fraction at 
z~lof^21± 6%, again with little variation across the 
luminosity and stellar mass bins. While the interloper 
fractions at both z ~ and z ~ 1 are consistent within 
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Fig. 6. — Top panels: Satellite velocity dispersion, measured within R p = [20, 150] kpc, for host galaxies at z ~ 1 (diamonds) and 
z ~ (circles) from sample A. In the left panel we plot the dispersion as a function of the S-band absolute magnitude of the host galaxy, 
while the right panel is as a function of stellar mass (M„). In addition we show the z ~ 1 measure men ts modified by a one magnitude 
dimming in the B-band luminosity of host galaxies (stars). The dashed lines are from Equations ITfl and llUI Bottom left panel: Mass-to-light 
ratios as a function of host galaxy Mg. Bottom right panel: Virial-to-stellar mass ratio as a function of host galaxy stellar mass. Since 
virial mass is defined with respect to the redshift-dependent critical density, an increase in M200/-M* of 1.3 is expected for an intrinsically 
non-evolving M20Q/-M* — M t relation (the same is true as a function of Mg). Symbol types are the same as the top panels. 



la, z ~ 1 host galaxies do have a slightly higher fraction 
of interlopers than z ~ hosts. A higher interloper rate 
at higher redshift might be attributable to the fact that 
more systems were still in the process of assembling then, 
and hence the higher fraction of interlopers could be re- 
flecting a higher fraction of systems that have not yet 
settled into dynamical equilibrium. Such a trend should 
not bias our velocity dispersion and mass estimates be- 
cause the interlopers are effectively accounted for in the 
dispersion measurements, regardless of their fraction. 15 
Figure includes power-law fits to the observed a — Lb 
and a — relations: 
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15 These interloper fractions are likely underestimated be- 
cause we have assumed that interlopers are distributed uniformly 
in dV while tests based on mock galaxy c atalogs suggest that 
the i n terloper distributi on is more complex Ivan den Bos ch et alj 
[2004]; |Chen_gtal] |20 03). Based on tests with mock catalogs, 
Ivan den Bosch et alj 120041) found that the assumption of a con- 
stant distribution of interlopers underestimates the interloper frac- 
tion by as much as 50%. Note, however, that these details do not 
impact the recovered satellite velocity dispersion and will affect the 
interloper fractions at z ~ 1 and z ~ in the same way. 



with la errors on the exponents of 0.1. We use only the 
lower three a(M*) points for the fit at z ~ because 
the dispersion in the highest stellar mass bin seems to 
deviate strongly from the relation indicated by the other 
points. Although the a — Lb relation at z ~ 1 is formally 
steeper than the corresponding relation at z ~ 0, we 
emphasize that the two slopes are consistent with one 
another, and stress that the inferred slope of the a — Lb 
relation depends on the range of luminosities probed and 
hence the exponents in this relation should in general be 
treated with caution (the same can be said of the a — M* 
relation) . 

For an NFW density profile with no anisotropy, virial 
mass is related to velocity dispersion as M200 oc a 2 ' 5 at 
Rp — 100 h^ 1 kpc (which is roughly the mean satel- 
lite R p used in our analysis), and thus M200/LB oc 
a 2 ' 5 J Lb- Hence, at z ~ 0, Equation implies that 
M 2 oo/L B oc £2.5x0.4-1 ^ constant, while at z ~ 1 
M 2 oo/L B oc L°- 5±0 - 3 . Notice that in general, if a oc L a , 
M/L will increase with L only when a > 0.4. Similar 
equations can be derived as a function of host galaxy 
stellar mass. These inferred M 2 oo/Lb versus Lb trends 
are perfectly consistent with the observed trends in the 
bottom panels of Figure HO In addition, the slope of our 
measured a — Lb relation at z ~ is c onsistent with 
previous estimates from satellite dynamics l|McKav et alJ 
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Fig. 7. — Satellite velocity dispersion as a function of host galaxy 
U — B color and Mb for samples at z ~ (circles) and z ~ 1 
(diamonds). At fixed Mg there is a clear difference between all 
(solid, black), red (open red), and blue (open with crosses, blue) 
hosts both at z ~ 1 and z ~ 0. 
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One should keep in mind that the er — Lb and cr — M* 
relations are not directly comparable because a bin in 
Lb contains a different fraction of red host galaxies from 
a similar bin in stellar mass. This is due to the fact that 
stellar mass is much more strongly correlated with galaxy 
color than i?-band absolute magnitude. At z ~ the 
fraction of host galaxies with red U — B color increases 
from 48% to 77% in our lowest to highest B-band mag- 
nitude bins, while the red fraction increases from 26% to 
91% from the lowest to highest stellar mass bins. Simi- 
larly, at z ~ 1 the fraction of red host galaxies increases 
from 33% to 55% from the faintest to brightest B-band 
magnitude bin, and from 4% to 88% for the smallest to 
largest stellar mass bins. The strongly varying fraction 
of red host galaxies as a function of host stellar mass 
makes the comparison and interpretation of the evolu- 
tion in <t(M*) between z ~ 1 and z ~ simpler because 
each end of the a — M* relation is dominated by a sin- 
gle population (red galaxies at the massive end and blue 
galaxies at the faint end). 

The a — M* and M200/M* — relations show little 
evolution between z ~ 1 and z ~ except for the high- 
est stellar mass bin. For host galaxies with M* < 10 11 
h~ 2 M Q the raw increase in M200/M* is constrained to 
1.4 ± 0.9 (but see below). For host galaxies with larger 
stellar mass the virial-to-stellar mass ratio increases by 
a factor of 4 ± 3 between these epochs. These hosts 
have large dispersions (> 300 km s _1 ) that correspond 
to massive dark matter halos (> 10 /i _1 M ), which 
commonly contain groups of galaxies, and hence they 
probably should not be interpreted in the same way as 
less massive host galaxies. We defer a more detailed dis- 
cussion of this extreme population to H7.2I 

When interpreting the evolution in the virial mass-to- 
light ratio, one should keep in mind that our definition 
of mass changes with redshift since the critical density is 
rcdshift-dependent. Specifically, at a fixed velocity dis- 
persion, the virial mass defined according to a region 
enclosing a mean density 200 times the critical density, 
M200, increases by a factor of ~ 1.3 from z ~ 1 to z ~ 0. 



A changing virial mass with redshift, even for a static, 
intrinsically non-evolving dark matter halo (i.e., a halo 
that is not accreting new material) , is a generic feature of 
all common virial definitions. In other words, an increase 
in M200/M* of 1.3 between z ~ 1 and z ~ is expected 
for an intrinsically non-evolving M200 /-M* — M* relation. 

Such an increase is precisely what is observed for host 
galaxies with M* < 10 11 h~ 2 M Q , which show a raw in- 
crease in Miqq/M* of 1.4 ± 0.9. Hence even this small 
increase does not reflect physically interesting changes 
in host galaxies. Folding in the changing virial defini- 
tion with redshift results in an "intrinsic" increase in 
AI200/M* of only 1.1 ± 0.7. This differentiation between 
intrinsic and definitional changes in M200/M* is sup- 
ported by the a — M* relation for host galaxies, which 
displays a very small raw increase from z ~ 1 to z ~ 
for hosts with M* < 10 11 h~ 2 M Q . Below, the intrinsic 
growth in M200/-W* with time will refer to growth after 
removing the factor of 1.3. 

To facilitate comparisons between z ~ 1 and z ~ 0, 
in the left hand panels of Figure we additionally plot 
the satellite velocity dispersion and M200/L for galaxies 
at z ~ 1 where all Mb values have been dimmed by 
one magnitude to account for evolution in M* between 
the medi an redshifts of the two samples (= 1.37 x (0.84 — 
0.06); cf. lFaber et al.ll2005D . This brings the two samples 
into near agreement. Our results as a function of Mb are 
thus consistent with host galaxies at z ~ 1 evolving into 
host galaxies at z ~ if their luminosities on average dim 
by one magnitude and their dark matter halo masses do 
not grow between these epochs. However, we emphasize 
that since this interpretation neglects the growth of dark 
matter halos with time it is rather unrealistic, and return 
to a more detailed interpretation of these results in H7.2I 

5.2. Results as a Function of Host Galaxy Color 

Finally, we investigate the a— Lb relation as a function 
of host galaxy U—B color. The a— M* relation cannot be 
probed as a function of color because the high/low stel- 
lar mass bins are almost entirely dominated by red/blue 
galaxies, so it is impossible to separate the host popu- 
lation into red and blue at these extremes of the stel- 
lar mass distribution and obtain adequate statistics, but 
these effects are more modest for a vs. Mb- Figure 
[3 and Table 01 present the satellite velocity dispersion 
as a function of host galaxy luminosity, color, and red- 
shift. We can only probe one broad bin in host luminosity 
(-20.5 < M B - 51og(ft) < -22.0) at z ~ 1 due to the 
limited number of satellites available. For every bin in 
host .B-band absolute magnitude there is a clear depen- 
dence of satellite velocity dispersion on host U — B color, 
at both z ~ 1 and z ~ 0. 

6. COMPARISON TO A SEMI-ANALYTIC MODEL 

We now compare our results to predictions from 
a semi-analytic model (SAM) of galaxy formation 
base d on the Millennium Run (MR) iV-body simula- 
ti on llSnr ingel et all I2005D . which is described in detail 
in lCroton et alJ l) 2006L In brief, the model evolves the 
contents of dark matter halos self-consistently accord- 
ing to a set of simple physical prescriptions derived from 
observational and theoretical phenomenology that gov- 
ern the evolution of baryons in a cosmological setting. 
These prescriptions track a wide range of physics, in- 



13 



eluding the cooling of baryons, star-formation in galac- 
tic disks and merger induced starbursts, supernova and 
AGN feedback, the tidal stripping of gas, and black hole 
growth. This model matches an array of observational 
results at z ~ 0, including the bj— and if— band galaxy 
luminosity functions and the color bimodality visible in 
the color-magnitude relation. 

In order to a facilitate direct comparison to our obser- 
vational results, mock surveys have been constructed out 
of the MR that match the geometry and sampling rates 
of both the DEEP2 and SDSS surveys. We then apply 
exactly the same search criteria and analysis methods to 
these mock surveys as to the data. 

Figure [S] presents a comparison of M200/M* between 
the MR SAM and data at z ~ (top panel) and z ~ 1 
(bottom panel) 16 . The Maoo/M* - M* relations for MR 
host galaxies have shapes quite similar to the observed 
relations at both low and high redshift (note that the 
highest stellar mass bin at z ~ 1 in the MR has only 13 
satellites and hence that measurement is unstable). In 
addition, the normalization of this relation in the MR 
is in agreement with observations at z ~ 1, though sys- 
tematically lower compared to observations at z ~ 0. 
This qualitative agreement highlights the power of cur- 
rent generation N-body simulations, including the MR, 
which explicitly follow sub-structure within each virial- 
ized dark matter halo. Such sub-structures host the pop- 
ulation of satellite galaxies from which our analysis is de- 
rived, and are clearly important for accurate modeling of 
the internal dynamics of group and cluster systems. 

The MR also provides an explicit test of our method- 
ology since we can compare the halo mass inferred 
from satellite kinematics to the true average host dark 
matter halo mass, which is readily available from the 
MR. The good agreement in Figure |S1 between the true 
and satellite-derived masses is very encouraging, though 
not unexpected based on previous tests of the satellite 
methodology (cf. Appendix Q. The agreement is less 
than ideal at the highest stellar mass bins at both z ~ 1 
and z ~ 0, though at z ~ 1 the discrepancy is of limited 
significance due to the modest number of pairs in the MR 
in that bin. 

At z ~ the discrepancy between true and satellite- 
derived halo masses is potentially more interesting. At 
this late epoch the MR SAM contains many satellites 
whose host subhalos have been stripped below the res- 
olution limit of the simulation (the so-called "orphan" 
population). In this situation the MR implements a 
prescription for the dynamical evolution of such galax- 
ies rather than tracking its evolution explicitly with the 
subhalo. When such galaxies are removed from the MR, 
the agreement between true and satellite-derived masses 
in the two highest stellar mass bins at z ~ significantly 
improves while simultaneously preserving the agreement 
at lower stellar masses and higher redshift. This issue 
clearly warrants further study. 

7. DISCUSSION 
7.1. Compatibility with Weak Lensing Results 

16 The MR galaxy formation model was run with a Salpeter IMF; 
in order to compare to the observational stellar masses, which were 
computed with a Chabrier IMF, the MR stellar masses have been 
lowered by the known offset between these IMFs, 0.3dex. 
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Fig. 8. — Observed and predicted virial-to-stellar mass ratios, 
M20o/M t , at z ~ (top panel), and z ~ 1 (bottom panel). 
The observed virial-to-stellar mass ratios (circles) are compared 
to both the true average virial-to-stellar mass ratio (triangles) in 
the MR semi-analytic model and the ratio derived from applying 
our satellite kinematics measurement techniques to that sample 
(diamonds). 



Our results at z ~ are most directly compara ble to 
the weak lensing results of iHoekstra et alJ 1)20051) who 
measured the lensing signal around a sample of isolated 
galaxies observed by the Red-Sequence Cluster Survey. 
They found that M v i r /L oc L 5 for isolated galaxies 
brighter than M B - 51og(/i) = -20.0, and M vil /L ~ 
constant for fainter galaxies. Our results in Figure HO 
are consistent with these findings. In particular, note 
that our derived M^qo/Lb appears to flatten out at 
Mb < —20.0, and we find a best-fit power law exponent 
of 0.5 for M200/LB versus L B for the three data points 
brighter than Mb = —20.0 (though taken together, our 
data are consis t ent w ith a slope of 0.0). As discussed in 
IHoekstra et al.l l)2005|) . their results are in broad agree- 
ment with other weak lensing studies. This overall con- 
sistency is quite encouraging since these two methods for 
deriving virial masses are almost entirely independent. 

In addition there is rough agreement in the normal- 
ization of the relation s derived from the two techniques. 
IHoekstra etafl l|2005[) found M vir = 1.9 x 10 12 h~ 1 M Q 
for galaxies with luminosity Mb — 5\og(h) = —20.6 at 
z = 0.32. In order to compare our results to theirs, 
we have evolved our absolute magnitudes by 1.37-2 to 
their median redshift, converted our definition of mass 
to theirs (their definition of virial mass is a region that 
encloses a mean density A V j r times the mean density of 
the universe, where A v ; r = 337 at z — and varies with 
redshift), and corrected for the ~ 30% over-prediction in 
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our derived mass due to incompleteness (cl. Appendix 
EJ). With these corrections we find a virial mass of 
M v ir = 1-8 x 10 12 Ii~ 1 Mq for host galaxies with a mean 
absolute magnitude Mb — 51og(/i) = —20.4. While this 
agreement is encouraging, at lower luminosities our cor- 
rected masses t end t o be higher than those quoted in 
iHoekstra et~aTI <|2005j) . by about 40% (though they still 
a gree within 2a) . 

IMandelbaum et all (|20Q6l M05) measured the aver- 
age virial-to-stellar mass ratio for galaxies in the SDSS 
as a function of galaxy stellar mass and found that 
M200/M* = 37 ± 13 h (47 ± 14 h) for early- type galax- 
ies with log(M*) = 10.5(10.7) /i~ 2 M (confidence levels 
quoted in M05 are 95%). In the present work we find 
M200/M* = 64± 2 2 ' h (62± 2 ^ h) for host galaxies with 
log(M*) = 10.6(10.8) h- 2 M Q , at z ~ (after correct- 
ing for the ~ 30% overestimate in our masses due to 
incompleteness; cf. Appendix lA"|l . Our results are hence 
quite consistent with M05 at the 95% confidence level. 
In addition, including the stellar mass of satellites results 
in a virial-to-stellar mass ratio of the entire system for 
(Af 2 ooA^*)tot = 52tg° h for host-satellites systems with 
log(M*, tot ) = 10.6 h- 2 M G . 

Perhaps most importantly, our results agree with M05 
in the highest stellar mass bin. For host galaxies with 
log(M„) = 11.1 h~ 2 M Q we find M 200 /M* = 256l^g h 
(again corrected for the ~ 30% overestimate in A/ 2 oo) 
while M05 found M 20 o/M* = 284±^ h for early- type 
galaxies with me an stellar mass l o gfM,,. ) = 11.3 /i _2 M . 

Most recently, iHevmans et alJ l|2006l H06) measured 
the weak lensing signal for galaxies in the COMBO- 17 
and GEMS surveys. These authors found that the virial- 
to-stellar mass ratio decreases from z ~ 0.8 to z ~ by 
at most a factor of 2.6 (la confidence) and are consistent 
with a constant value of M 2 oo/-M* = 661 2 q h over this 
redshift range. These authors include only galaxies with 
log(M*) > 10.2 h~ 2 MQ, where there data is complete 
for z < 0.8. The evolution of M 2 oo/M* measured in 
the present work is again in good agreement with these 
weak lensing results, including the non-evolving virial-to- 
stellar mass ratio from z ~ 1 to z ~ for galaxies with 
M* < 10 11 Ii~ 2 Mq (the average stellar mass of galaxies 
in H06 is (log(M»)) = 10.5 h~ 2 M Q ). 

A more detailed comparison between these two meth- 
ods for estimating halo masses is complicated because 
the systems that satellite dynamics probe are systemat- 
ically different from the systems weak-lensing analyses 
explore. In order to uncover any measurement bias one 
would need to measure the weak-lensing signal for the 
same set of galaxies that are used in satellite dynamics 
studies, which is likely feasible given the current number 
of host galaxies found in the SDSS. 

7.2. The Evolution of Host Galaxies from z ~ 1 to z ~ 

We now explore the evolution in the virial mass-to- 
light and virial-to-stellar mass ratios i.e., A/ 2 oo/£b and 
Af 2 oo/M», respectively - between z ~ 1 and z ~ 0. Evo- 
lution in M 2 oo/M* is somewhat simpler to interpret be- 
cause stellar mass can only increase with time, and so we 
consider it first. 

Our data suggest a bimodality in the evolutionary his- 
tory of host galaxies between z ~ 1 and z ~ 0: host 
galaxies below M* ~ 10 11 h^ 2 M Q maintain a roughly 



constant virial-to-stellar mass ratio (the intrinsic ratio 
increases by a factor of 1.1 ± 0.7) while hosts above this 
stellar mass scale, which are predominantly red in U — B 
color, experience a factor of 4 ± 3 increase in M 2 oo/M». 
The quenching of star formation, whic h becomes particu- 
larly effective above this stellar mass llDekel fc Birnboiml 
l2006HCroton et alJl2006HCattaneo et alJl2006[) . is a kev 
physical process likely responsible for this bimodality. 

Galaxies below M» ~ 10 11 h~ 2 M Q continue to form 
stars between z ~ 1 and z ~ 0. Recent modeling of 
the star formation history of relatively low-mass blue 
galaxies has suggested that these objects grow in stel- 
lar mass by roughly a factor of 2 between these epochs 
ijNoeske et al.ll200a() . Cosmological simulations of dark 
matter reveal that ~ 10 12 h~ 1 M Q dark matter halos - 
such as the ones which our typical host galaxies appear 
to reside in - gro w on average by a fa ctor of ~ 2 over 
this time period ijWechsler et alJl2"002j) . The growth in 
stellar mass of host galaxies hence proceeds in lockstep 
with the growth of their dark matter halos, yielding a 
constant virial-to-stellar mass ratio from z ~ 1 to z ~ 0. 

If the virial-to-stellar mass ratio for these galaxies does 
not change with z, their star formation efficiency must 
also remain constant over this time. Specifically, after 
correcting for the ~ 30% overestimate in our masses 
due to incompleteness (cf. Appendix lA")l . our results im- 
ply a star formation efficiency for host galaxies of 77 = 

m, sip. o.l3, assuming a universal baryon fraction of 

M200 sjl ' J 

0.17 l|Snergel et alJ2006t) and h = 0.7. Including the stel- 
lar mass of the satellite galaxies yields an efficiency for 
the entire system of r\ as 0.16. These efficiency factors are 
bracketed by values of r\ = 0.08 derived from the g lobal 
baryonic mass function l|Read fc Trenthaml I2005D and 
various weak-lensing studies, which find ri ~ .05 — 0.3 
l|Mandelbaum et alJl2006t IHevmans et al.ll2006[) . 

The lack of evolution in M200/M* found here is 
broadly consistent with measurements of the evolution 
of the stellar mass Tully-Fishe r relation (TFR) over 
this time period. In partic ular, iConselice et alJ (2005) 
and lBoehm fc Zieglerl {2006) converted their stellar mass 
TFR into a dynamical-to-stellar mass ratio and found lit- 
tle evolution from z ~ 1 to z ~ 0; using ~ 550 galaxies, 
iKassin et al.l l)2006a[) found no evolution in the stellar 
mass TFR between these epochs, which under simple as- 
sumptions translates into a non-evolving dynamical-to- 
stellar mass ratio. One should keep in mind that there 
are numerous details which make a one-to-one compari- 
son between the stellar mass TFR and our results compli- 
cated. For example, the relation between internal galaxy 
kinematics (e.g. the maximum circular velocity of the 
disk) and dark matter halo kinematics is non-trivial, and 
likely changes with time. Recent SPH simulations of indi- 
vidual disk galaxies have found that stellar mass and cir- 
cular velocity co-evolve (both roughly doubling between 
z ~ 1 and z ~ 0), yielding a non-evolv ing stellar mass 
TFR (|Portinari fc Sommer-Larsenll2006ft . 

The evolutionary history of the most massive (M* > 
10 11 h~ 2 MQ, and therefore red) host galaxies is pre- 
dicted to be qualitatively different. A variety of recent 
models suggest that feedback mechanisms in high mass 
halos prevent the ambient gas from cooling and form- 
ing new stars |C attane o et al J l2*006i iDekel fc Birnboinl 
l2006HCroton et al.l2006fl . Hence the stellar mass of these 
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galaxies can only increase via merging with other galax- 
ies. Numerous observational results suggest that the stel- 
lar ma ss in bright red galaxies at most doubles since 
z ~ 1 (|Borch et al.112009: iBundv et al.ll2005al iBell et all 
2004). Yet massive dark mat ter halos are growing r apidly 
between z ~ 1 and z ~ ( Wechsl er et alJ 12002'). The 
quenching of star formation within these dark matter 
halos provides a natural, qualitative, explanation for the 
observed growth in M 2 oo/M* with time. 

The quenching of star formation in these massive sys- 
tems implies that their star formation efficiency decreases 
with time since new gas that comes into these halos is 
prevented from condensing to form new stars. The ob- 
served increase in M200/M* between z ~ 1 and z ~ 
for these systems corresponds to a decrease in efficiency 
from 14% to 3% between these epochs. 

This simple picture is complicated by our selection cri- 
teria. In particular, host galaxies are required to not have 
any comparably bright companions within a search cylin- 
der. This non-trivial selection makes it difficult to reason 
generally about the growth of halos and stellar mass with 
time, and favors comparisons to more detailed models of 
galaxy formation (see ^and below). However, these se- 
lection effects may also account for part of the observed 
increase in M200/-M* between z ~ 1 and z ^ 0. 

Due to the continual merging of galaxies within ha- 
los over time, galaxies can enter the host galaxy sam- 
ple between z ~ 1 and z ~ 0. In particular, objects 
within groups and clusters of galaxies at z ~ 1, which 
are in halos that are more massive than the average host 
galaxy dark matter halo, may not be classified as host 
galaxies due to their bright companions, but can turn 
into what we would call host galaxies at late times. As 
groups and clusters at z ~ 1 evolve, dynamical friction 
will cause some fraction of bright companion galaxies to 
merge with the central galaxy. Those systems in which 
only one luminous galaxy remains within a massive dark 
matter halo will be simil ar to so-called "fossil groups" 
(see e.g. Uones et al]l2003|) . The increase of fossil groups 
between z ~ 1 and z ~ is a possible explanation for 
the increase in M200/-W* for high-stellar-mass host galax- 
ies. This mechanism will likely not affect the evolution 
of M200/M* at smaller host stellar masses because the 
merging of ~ L* galaxi es in massive halos will gener- 
ate massive red galaxies (fHopkins et all2Q 06) , not lower- 
mass, blue galaxies. 

The stellar mass and redshift dependence of the virial- 
to-stellar mass ratio of host galaxies is reproduced qual- 
itatively in the M illennium Run semi-analytic model 
l)Croton et aU feOOfi). In particular, this model manages 
to capture the stellar mass dependence of M200/M* at 
z ~ 1 quite well. At z ~ this model reproduces the ob- 
served shape of the M 2 oo/M* — M„ relation, though with 
a lower normalization. The Millennium Run, and mod- 
els like it, are particularly useful for interpreting these 
results because these semi-analytic models track the evo- 
lutionary history of mock galaxies. Hence, for example, 
we are capable of directly confronting the explanations 
presented above with a self-consistent model. We can ask 
where host galaxies at z ~ 1 end up at z ~ and, con- 
versely, where z ~ host galaxies were at z ~ 1. Such a 
comparison will be the focus of future work. 

Evolution in the relation between host galaxy virial 
mass and £?-band absolute magnitude is more difficult to 



interpret because Mb can both increase, due to recent 
episodes of star-formation and mergers, and decrease, 
due to fading in the stellar population. Although it was 
demonstrated earlier that our results as a function of Mb 
are consistent with the average host galaxy Mb fading 
by one magnitude from z ~ 1 and z ~ while its dark 
matter halo does not evolve, our results as a function 
of stellar mass cast this simple interpretation into doubt. 
In particular, evolution in the M200/-M* relation strongly 
suggests that the average halo mass of host galaxies is in- 
creasing from z ~ 1 to z ~ 0, since the stellar mass of 
host galaxies is most likely increasing with time (and, 
moreover, all currently favored cosmological models in- 
dicate that halos grow in mass with time). If halo mass 
increases significantly, host galaxies at z ~ 1 must fade 
by less than one magnitude (i.e. fading by less than 
L* between these epochs) in order to reach the z ~ 
M200/LB relation (cf. Figure HJ left panels). 

8. SUMMARY 
We summarize our main results and conclusions: 

1 . The U — B color distribution of satellite host galax- 
ies is indistinguishable from the colors of all galax- 
ies (of similar luminosities), at both z ~ 1 and 
z ~ 0. Satellites of host galaxies at z ~ 1 are 
on average slightly fainter in the -B-band relative 
to their host galaxy compared to z ~ 0, and host 
galaxies on average have a comparable number of 
satellites at both epochs. 

2. The line-of-sight velocity dispersion, a, of satel- 
lites increases with the average host galaxy lumi- 
nosity, Lb, and stellar mass, M*. In particular, 
cr oc M°- 4±al at both z ~ 1 and z ~ 0, while 
cr cx L° B 6±01 at z ~ 1 and a oc L B A±0A at z ~ 
(though it should be noted that these exponents 
are sensitive to the precise range of Lb and M* 
over which the relation is fit). In addition, at fixed 
Mb red host galaxies have larger satellite velocity 
dispersions compared to blue hosts, both at z ~ 1 
and z ~ 0. 

3. The virial-to-stellar mass ratio, M200/M*, for host 
galaxies with M* < 10 11 h^ 2 M & increases by a 
factor of 1.1 ± 0.7 from z ~ 1 to z ~ 0; the ratio 
does not evolve significantly between these epochs. 
Host galaxies with stellar mass above M* ~ 10 11 
h Mq, which are predominantly red in restframe 
U — B color, experience an increase in M200/M* of 
a factor of 4 ± 3 between these epochs. 

4. The Millennium Run semi-analytic model of 
galaxy evolution reproduces the observed trends of 
M200/M1, as a function of M* at z ~ 1 quite well, 
while at z ~ the model reproduces the shape 
of the M200/-M* — M* but with a somewhat lower 
normalization compared to observations. 
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APPENDIX 

A. IMPACT OF SURVEY COMPLETENESS ON THE DERIVED AVERAGE HOST GALAXY HALO 

MASS 

As described in H3.2I the DEEP2 survey obtains redshifts for ~ 50% of all galaxies with R < 24.1. In an incomplete 
galaxy catalog, we may identify a galaxy as isolated when in reality it has a bright companion that simply failed to 
be targeted. We assess the impact of incompleteness on the recovered velocity dispersion and halo mass by randomly 
diluting the SDSS survey and then searching for host galaxies and satellites. We show in Figure lA"9l the change in the 
measured halo mass, M200, as a function of Mb and M*, comparing host and satellite galaxies identified in the complete 
SDSS to the diluted SDSS. The recovered masses agree within ~ lcr, which is encouraging, as this demonstrates that 
the modest incompleteness of the DEEP2 survey does not greatly change the recovered dispersion. 

There is, however, a noticeable bias in the recovered masses: the diluted samples yield higher masses at a given host 
galaxy Mb and M*. This is due to the inclusion of systems that are not truly isolated, but in fact reside in denser 
environments, which will have on average larger masses than isolated systems. We find that ~ 35% of host galaxies 
identifie d in the diluted sampl e are in fact not truly isolated. This number is comparable to the contamination rate 
found in iConrov et al J l)2005b|) . who analyzed mock catalogs that included the same selection effects as the DEEP2 
survey. This contamination fraction depends mildly on host luminosity, with sub-samples of brighter hosts having 
a lower contamination fraction than sub-samples of fainter hosts. This trend reflects the fact that brighter galaxies 
tend to be found in denser environments, so they are more likely to still have bright companions after dilution. The 
numerous bright companions will help prevent a bright galaxy from being falsely identified as isolated in the diluted 
sample, while the relative paucity of similar luminosity companions around fainter galaxies will make the identification 
of such a galaxy as (falsely) isolated in a diluted sample easier. 

The aver age ratio between diluted and complete SDSS host galaxy halo masses is indicted by the dashed lines in 
Figure EH We use these average mass corrections to arrive at host galaxy halo masses for 100% complete surveys, 
both at z ~ and z ~ 1, when comparing against other work. 

These falsely isolated host galaxies probably do not strongly effect the recovered velocity dispersion because the 
satellites that they contribute to the dV distribution resemble the interloper population: they are a small fraction 
of the total sample and they are more uniformly distributed in dV compared to satellites from truly isolated host 
galaxies. 

Finally, we point out that any potential bias that may result from incompleteness will only impact the absolute a 
and M200/L values quoted herein. Results concerning the relative change in these quantities between z ~ 1 and z ~ 
are likely not affected by incompleteness, since we have diluted the SDSS down to the same completeness as DEEP2 
for those studies. 

B. IMPACT OF CHANGING MEASUREMENT PARAMETERS 

In this section we explore the effects of varying the parameters we have used when measuring velocity dispersions 
and deriving virial masses. Specifically, we explore the impact of changing the host-satellite search criteria or the 
number and size of the radial bins in which we measure velocity dispersions, and in addition we show that the derived 
masses are insensitive to the assumed halo concentration and velocity anisotropy. For sim plicit y, in the following tests 
we have selected one sample at z ~ and one at z ~ 1 to illustrate our results (see Table lB5|l . and the z ~ sample 
has been diluted to match the completeness of the sample at z ~ 1, as mentioned previously. Our fiducial set of 
assumptions and parameters are an isotropic velocity distribution (/3 — 0), an NFW concentration of c = 10, a bin in 
projected separation between satellite and host of R p = [20, 150] hr 1 kpc, and search criteria A. 

We begin by te sting assumptions that affect the measurement of velocity dispersions (and through this, the derived 
masses). In Table lB5l we show the change in the derived virial mass when wc use different search criteria and different 
radial bins. The primary effect of changing search criteria is that the number of satellites changes, which in turn alters 
the errors on the recovered velocity dispersion and virial mass. For example, using search criteria B or C at z ~ 1 
result s in only 30 and 32 satellites within R p — [20, 150] /i -1 kpc, respectively, for the luminosity bins we consider in 
Table lB5l hence the velocity dispersions and masses estimated for these samples are rather unstable. As can be gleaned 
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Fig. A9. — Ratio of host galaxy halo masses derived from the complete SDSS survey to masses derived from the SDSS survey diluted by 
40% to match the completeness of the DEEP2 survey. The diluted sample yields systematically higher virial masses as a function of both 
Mg (left panel) and M» (right panel). The dashed line in each panel represents the average 
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for details. 



TABLE B5 

Dependence of derived virial mass on various parameters 



Sample 


Anisotropy 
/3 = ft 


Concentration [ 
c = 5 c = 10 c=15 


A 


Search Criteria [ 
B C 


Radial bins (h 
[20, 100] [20, 150] 


1 kpc) 
[20, 200] 


SDSS 
DEEP2 


2.8 ±0.6 
1.5 ± 0.6 


2.8 ±0.6 
1.6 ±0.7 


2.8 ±0.7 2.8 ±0.6 2.8 ± 0.6 
1.3 ±0.6 1.5 ±0.6 1.6 ±0.6 


2.8 ±0 
1.5 ± 


6 3.0 ±0.8 3.0 ±1.0 
6 1.3 ±1.0 2.5 ±2.5 


2.3 ±0.6 2.8 ±0.6 
1.0 ±0.5 1.5 ±0.6 


3.5 ±0.7 
2.5 ± 1.0 



Note. — The SDSS and DEEP2 samples include host galaxies with -20.25 < M B - 51og(ft) < -19.75 and -21.0 < M B - 51og(ft) < 
— 19.5, respectively. All values are virial masses (M200) in units of 10 12 1 Mq. The fiducial set of assumptions and parameters are 
/3 = 0, c = 10, search criteria A, and radial bin [20, 150] h _1 kpc. Each table entry is the result of varying only one of these parameters 
at a time. See Tabled for details regarding the search criteria. 



from the table, changing the search criteria or varying the radial binning results in changes in the recovered masses 
within ler of our fiducial estimates. Note also that using search criteria B at z ~ 1 and C at z ~ amounts to using 
consistent isolation criteria in comoving coordinates (because in this case the projected separation scales as (l + z)" 1 ), 
and hence we conclude that using such coordinates, as opposed to physical coordinates, results in no significant change 
in our results. 

We now test the impact on a of using relatively bright satellites, i.e satellites that are only one magnitude fainter 
than their host galaxy. One might expect more massive (brighter) satellites to be out of virial equilibrium with the 
host galaxy since dynamical friction acts more efficiently on these satellites compared to less massive satellites. Using 
the complete SDSS sample, we measure the satellite velocity dispersion for hosts identified using search criteria A as 
a function of dM, the separation in brightness between the satellite and host galaxies. Specifically, we only include 
satellites that are at least 1.0, 1.5, 2.0, 2.5, and 3.0 magnitudes fainter than their host galaxy, and find that the 
dispersions estimated from these restricted samples all agree well within la. 

Since using different radial bins causes the largest changes in these tests, we have explored its effects in more detail. 
We have reconstructed the a — L and er — M* relations using each set of radial bins and find no qualitative change in 
our results. In particular the exponents of these relations when using R p = [20,200] h~ x kpc are the same at both 
z ~ 1 and z ~ 0, in agreement with the results shown in Figure El The results for the R p = [20, 100] h^ 1 kpc are 
noisier, as there are ~ 50% fewer satellites in this sample, but are also consistent with our fiducial results. 

We now test assumptions that only affect the conversion between velocity dispersion and mass. In Table lB5l we show 
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the effects of varying the concentration, c, or velocity anisotropy, f3, on the recovered virial mass. We use concentrations 
of 5, 10, and 15, which span the range of typical concentrations for ~ 10 12 Ii^Mq dark matter hal os for < z < 1. 
We te st two forms of the velocity anisotropy, the isotropic case (/3 = 0), and a form suggested bv iMamon fc Lokasl 
1)2005(1 that is a good fit to a compilation of anisotropies derived from clusters within iV-body simulations: 

Mr) = , (Bl) 

2r + r a 

where r a — 0.18 r2oo- This functional form has not been tested in galaxy-sized halos, though see the discussion in 
Appendix [Q 

At first glance it seems somewhat remarkable that the derived virial mass is almost independent of the concentration 
and velocity anisotropy parameters over the range tested. As it turns out, the line-of-sight velocity dispersion profile in 
the regime 50 < R p < 150 hr 1 kpc, which encompasses most of the satellites we use, is almost completely insensitive 
to concentration and anisotropy, while on both smaller and larger scales their effects are quite noticeable. We are, by 
necessity if not luck, probing a "sweet spot" in the velocity dispersion profile. 

In addition, we have tested the sensitivity of the derived virial mass to the number of radial bins used to measure 
velocity dispersions. We can only perform this test for the SDSS sample due to the limited size of the DEEP2 sample. 
We have measured the velocity dispersion in three bins of projected separation, R p = [20 — 100], [100 — 200], and 
[200 - 300] h- 1 kpc, for the SDSS sample used in Table lB5l We find a best fit virial mass of M 200 = 3.8 ± 0.5 x 10 12 
fe _1 M0, which is within la of the values quoted in Table lB5l using only one radial bin. Furthermore, rebinning the 
undiluted SDSS sample results in a mass of Af200 = 2.6 ± 0.6 x 10 12 fe" 1 ]^ " ^. Wh ile we see a decline in the velocity 
dispersion profile for the undiluted sample in agreement with lPrada et all l(2003fl . we see no decline for the diluted 
sample. This is likely due to the fact that the 35% of host galaxies in the diluted sample that are not truly isolated 
introduce noise into the dispersion measurement on large scales. 

Finally, we turn to our assumption of an NFW l|Navarro et alJll"997() characterization of the density profile. Recently, 
iPrada et alJ l(2003f l have detected a decline in the satellite velocity dispersion profile with increasing distance from the 
host galaxy at z ~ 0. This decline distinctly favors an NFW density profile over an isothermal distribution. However, 
we cannot distinguish between these and other density distributions at z ~ 1 because there are an insufficient number 
of satellites, making it difficult to bin more finely in radius and probe the velocity dispersion profile. 

We are motivated to choose the same parameterization of the density profile at z ~ 1 and z ~ because it 
seems unphysical for the density profile to drastically change, e.g. from isothermal to NFW, over this interval. 
Furthermore, an NFW-like density profile (notwithstanding minor deviations) appears to be a gen eric feature of 
hierarchica l clustering, and is not sensitive to initial conditions or cosmological parameters (see e.g. iNavarro et alJ 
Il997il200l . Note that our conclusions regarding the dependence of the satellite velocity dispersion on host luminosity 
and redshift are completely independent of density profile considerations. 

C. ASSESSMENT OF ASSUMPTIONS 

We identify three main assumptions inherent in using satellite galaxies to probe the virial dark matter mass of their 
hosts: 1) there is little scatter between host galaxy luminosity and dark matter halo virial mass, 2) satellites are 
fair tracers of the underlying dark matter velocity field, and 3) the velocity difference distribution of satellites and 
interlopers can be modeled as a Gaussian and a constant, respectively. These assumptions are well motivated both 
observationally and theoretically. 

The first assumption (low scatter between mass and luminosity) must hold for our stacking procedure to work; 
it is motivated by the Tully-Fisher relation for disk-dominated galaxies and the Faber- Jackson relation for bulge- 
dominated galaxies. These relations have only a modest amount of scatter; typically 0.2 — 0.9 magnitudes at fixed 
velocity. The scatter is even sma ller when these relations are quoted as a function of stellar mass or total baryonic mass 
l|Verheiienll2001i fMcGau gh 2005). These tight correlations, with the assumption that the velocity measured probes the 
underlying halo mass, suggest that the scatter between halo mass and galaxy luminosity is not large, and hence that 
our stacking procedure is valid. In addition, models with a ti ght relation between galaxy luminosity and halo mass 
successfully reproduce a wide variety of observa tions (see e.g. iCole et al.ll2000t iTasitsiomi et aljf2004t iConrov et alJ 
I2005at ICroton et al]l200fil IVale fc Ostrikei1l200l . 

Under the assumption that satellite galaxies can be identified with subhalos in dissipationlcss A^-body simulations, 
these simulations support our second assumption, that satellites are fair tracers of the underlying mass distribution. 
A number of studies have looked for a "velocity bias" between dark matter and subhalos. Recent work has suggested 
that such a bias, defined as the r atio between the velocity dispersion of s ubhalos and d ark matter, is about 10% 
when averaged over entire clusters <)Ohigna et al]l200(l iDiemand et al.l l2004'). In fact, when iFaltenbacher fc Diemandl 
( 2006) identify subhalos in a simulated cluster in such a way that their spatial distribution matches observed galaxies, 
they find no veloc ity bias. Although no study has systematically investigated the velocity bias in galaxy-size halos, 
IPrada et alJ l(2003j) showed that subhalos identified in a simulated galaxy-sized halo accurately reflect the underlying 
mass distribution. Admittedly, the velocity bias of the most massive subhalos, which should correspond to the satellites 
studied here (since we use relatively bright satellites, which are likely associated with the most massive subhalos), has 
not been adequately investigated. Nevertheless, we conclude that satellite galaxies, should be fair tracers of the 
underlying dark matter halo mass to ~ 10% or better. 

The second and third assumptions have been tested in conjunction by a number of authors. IPrada et alJ l|2003f) 
identified satellite galaxies with subhalos in dissipationless A^-body simulations, and artificially added in interloper 
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galaxies uniformly in phase space. They found that for an isolated Milky- Way sized dark matter halo, the Gaussian 
plus constant parameterization accurately recovers the velocity dispersion profile of subhalos, which in turn accurately 
r eflects the underlying dark m atter halo mass. 

Ivan den Bosch et alJ (2004) also tested these assumptions using cosmological mock galaxy catalogs at z ~ 0, and 
found that the simple Gaussian plus constant parameterization accurately recovers the velocity dispersion of host 
galaxies as a function of galaxy luminosity. These authors additionally investigated the impact of orbital anisotropy 
and a spatial (anti-) bias of the satellite galaxies, and found the resulting effects on the recovered dispersion to be 
minimal. Specifically, when using satellites within one-third of the virial radius from the host galaxy, they found that 
if j3 changes from —0.5 to 0.5 then the velocity dispersi on changes by < 10% . 

These assumptions have also been tested at z ~ 1 bv lConrov et alJ <l2005bt) who used cosmo logical simulations into 
which mock galaxies were inserted using a halo model approach (see lYan et alJ 120031 12004. for details concerning 
these mock catalogs). These authors found that the velocity dispersion profile of mock galaxies accurately reflects 
the underlying halo mass, and that interlopers can be effectively modeled as a constant component to the velocity 
difference distribution. 

Consistency between weak-lensing and satellite measurements of M200 1 L in past studies gives further confirmation 
that satellites are fair tracers of the underlying d ark matter velocity field (unless both satellite dynamics and weak 
lensing measurements are biased in a similar way). iMcKav et alJ (|2002) compared the dependence of the mass-to-light 
ratio on luminosity derived from both satellite dynamics and weak lensing measurements, and found agreement within 
Iff. Furthermore, both weak l ensing studies l|Guzik fc S gJjal^JgJjQglJH 

120051 iMandelbaum et al l 120061) and satellite dynamics i|McKav et al.l |2002t iPrada et al.l 120031 iBrainerd fc Snecian 
2003; van den Bosch et al.ll2004j) have found that the derived virial mass scales with host luminosity as M oc L a with 
1-0 S a S 1-5- The range in a can be at least partially explained by the different regimes (isolated galaxies vs. groups 
and clusters) and photometric bands probed. 



